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Klein paradox for arbitrary spatio-temporal scalar potential and Josephson-like 

current in graphene 
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We derive the exact time evolution according to the Dirac-Weyl equation, describing a mono- 
layer of graphene, in the presence of a scalar potential U{x, t) of arbitrary spatial and temporal 
dependence at normal incidence, Py = 0. This solution shows that the Klein paradox (the absence 
of backscattering) persists even for arbitrary temporal modulations of the barrier. Moreover, we 
identify an unusual oscillating current jy running along the barrier, despite of the vanishing momen- 
tum in y-direction. This current exhibits resemblence to the Josephson current in superconductors, 
including the occurance of Shapiro steps and its sine-like dependence on the phase difference of wave 
functions. 
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The growing interest to graphene, see e.g. Ref. Ill, |2;, is 
stimulated by many unusual and sometimes counterintu- 
itive properties of this two dimensional material. Indeed, 
graphene supplies charge carriers exhibiting the pseudo- 
relativistic dynamics of massless Dirac fermions. One 
example of the unusual dynamics of electrons and holes 
in graphene is the Klein tunneling phenomenon [3| which 
occurs with unit probability through arbitrarily high and 
thick barriers at perpendicular incidence, irrespective of 
the particle energy, in accordance with experiment [J|. 
In consequence, the question arose of how to control the 
electron motion in graphene and hence boosted detailed 
studies of Dirac fermions under the influence of various 
forms of scalar |5l4l2| or vector [l3| potentials. 

So far, lots of works were devoted to studies of 
graphene subject to static periodic electric fields, since 
these structures known as graphene superlattices p, |6[ 
allow controlling both specrum and transport properties 
of electrons in graphene. For instance, it was shown 
[a| that ID graphene superlattices have a deep anal- 
ogy with photonic crystals formed by alternating right- 
handed and left-handed transparent media, similar as the 
earlier stated analogy [7| of a p-n junction in graphene 
to a Veselago lens. Superlattices of electrostatic periodic 
potentials can be used to coUimate the directional spread 
of electron beams in graphene [8| so that waves of zero 
transverse momentum will dominate transport proper- 
ties. 

Applying a time-dependent laser field to a pristine 
graphene sample opens an alternative and efficient way 
[9l4ll| to control spectrum and transport properties of 
graphene samples. It has been shown that changing the 
time dependence of laser fields can mimic [10| the influ- 
ence of any electrostatic graphene superlattices on the 
electron spectrum in graphene. Studies of how electron 
transport in graphene is affected by time- and- space de- 
pendent potentials are yet limited. 

For example, it was shown [lO| that Dirac fermions in 
graphene superlattices can aquire an effective mass pro- 



portional to the frequency of an applied laser field, ac- 
companied with an exponential suppression of chiral tun- 
neling even for perpendicular incidence upon the barrier, 
similar as an ac-electric field directed parallel to the bar- 
rier pjj , in stark contrast to Klein tunneling occuring in 
the absence of the laser field. The question arises whether 
time-dependent modulations of a electrostatic barrier it- 
self would generate backscattering or not. In this Letter 
we shall answer this question and further demonstrate 
novel transport phenomena occurring for general space- 
time dependent potentials. We focus on a geometry ex- 
hibiting the most pronounced Klein paradox, at vanish- 
ing transverse momentum {py = 0) parallel to a barrier, 
taking the barrier as homogeneous along the y-direction. 
We present the general solution of the Dirac equation 
for a scalar potential U{x, t) of any spatial and temporal 
dependencies. For this case we demonstrate unit trans- 
mission probability, thus, generalizing the Klein paradox 
for arbitrary time dependence in U{x,t). Secondly, we 
discover an oscillating current jy along the barrier, which 
in view oi Py = seems counter intuitive. This current 
bears resemblance to the Josephson current in supercon- 
ductors and we demonstrate even similar properties, such 
as Shapiro steps for properly chosen frequencies. Exper- 
imental test of our predictions should be within reach of 
present day nanostructure design on graphene [12|. 

Model. — The honeycomb lattice of graphene engen- 
ders two copies, Tz — ±1, of Dirac-Weyl Hamiltonians 

H 



Hn = VF [r^a^Px + cTyPy] 



(1) 



centered about two inequivalent Dirac points ( "valleys" ) 
K and K' at corners of the hexagonal first Brillouin zone 
where electron-hole symmetric bands touch; the Pauli 
matrices cr^.y.z act on two-component spinors represent- 
ing sublattice amplitudes, exhibiting opposite Fermion 
helicities, -p/p = ±1 . SU(2) rotations with respect to 
the vector r of three Pauli matrices allow to continu- 
ously transform both copies into one another |15| which 



motivated the terminus "valleytronics" for isospin ma- 
nipulations based on eigenstates to Tz, Ref. 
ogy to the weh-known research area of spintronics 



161 . in anal- 
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Proposals exist to valley polarize carriers, by rneans of 
nanoribbons terminated by zig-zag edges 16J, [iSl, Il9| , by 
exploiting trigonal warping at elevated energies [20|, or 
by absorbing magnetic textures [21 1. 

Below, we focus on valley polarized situations. In- 
deed, smooth electromagnetic or disorder potentials do 
not couple the two valleys [22| so that calculations can 
be done independently, for either r^ = -1-1 or r^ = — 1. 
Including now the barrier potential U{x,t) = hW{x,t) 
the Dirac equation becomes 



v^iT^p^ + ipy)'i'A + hW{x,t)^'B =ifi- 



dt 
d-^B 



dt ' 

where wave functions ^^,^5 describe electrons on ei- 
ther of the triangular graphene sublattices, v-p is the 
Fermi velocity, and the momentum operator is defined 
as {px,Py) — {—ihd/dx,—ihd/dy). This equation has 
been solved analytically for time-independent potentials 
either by matching [3| of wave functions for rectangular 



barriers, or by the WKB method [23|-|25| for smooth bar- 



riers. Additional time dependent harmonic oscillations 
have been considered of gate voltages on either side of 
the rectangular barrier [26| , or of an electric field par allel 
to the barrier [ll| or in resonance approximation [lOj . In 
the two latter cases a dynamical gap opens which can 
suppress even perpendicular (py = 0) tunneling. 

Our goal here is to obtain the exact solution of eq. 
dD) for Py — and arbitrary potential W{x,t) acting 
at positive times, i.e. W{x,t < 0) = . In this case, 
the wave functions \1/a and ^b depend on time t and 
the coordinate x across the barrier, but not on the y- 
coordinate. This simplifies ^ to read 



Exact solution. — To solve (jSj we use the Ansatz 

1 / QiS±{x,t) 

^±(^'')=2Ur.e'^±(-*) 



(3) 



(4) 



where the sign ± distinguishes left and right propagating 
solutions. Inserting ^ into ^ results in 



dtS± {x, t) ± VpdxS± (x, t) + W(x,t)=Q 



(5) 



This first order partial differentia l eq uation can be solved 
by the method of characteristics [27[ , according to which 
the vector (1, ±wf, —W) is tangent to the solution S±, 



yielding 

S±{x,t) = s'i\xTVpt,Q)- / dt' W{xTVp{t-t'),t') 

Jo 

(6) 
explicitly in terms of W{x,t). In view of ([4]) the term 

S'^ (a;,0) describes the initial wave function ■0±(a;,O) at 
time t = which, in the absence of the barrier at i < 0, 
can be, e.g., a plane wave of wave number k in the x- 
direction, Sj. {x, 0) = ±kx, or some wave packets. Then 
eq. (O, together with (|4]), describes the full solution for 
t> 

^P{x,t) = a+{x - vpt) (^ \ e-'I>'' ^(— (*-*'),*') 

where a±{x) = e'"^± ^^^ = ['i'A{x,t = 0) ± T:,-^B{x,t = 
0)]/2 encodes the initial condition. In particular, if the 
wave packet is initially purely right moving, a_ = , 
according to ([7]), it continues propagating to the right 
at times t > with undistorted density distribution 
|a_|_(a; — v-pt)\'^ without reflection, acquiring at most a 
phase factor due to the time dependent barrier. How- 
ever, the situation becomes much more intriguing when 
we consider a superposition of left and right movers. 

Josephson-like current. — Now we evaluate the current 
density 



jx{x,t) = Vpip* {x,t)T^axip{x,t) 
3y{x,t) = Vj,'ip*{x,t)ayip{x,t) 



(8) 



for Py — 0. A very nontrivial result occurs if we consider 
a linear superposition 



ip{x, t) = a+%jj^{x,t) + a_^_(a;, t) 



(9) 



of the right and left moving solutions, -0+ and V'-j ob- 
tained above, eqs. (|4]) and ^, assuming that \a+\^ + 
|q;_P = 1/2 . Note, that an initial density peak arising 
from some voltage pulse will generally contain simulta- 
neously left and right moving amplitudes. 

For the x-component of the current density, flowing 
perpendicular to the barrier, we find 

jx{x,t) = Y {\a+a+{x - Vpt)\'^ - \a-a_{x + Vpt)\^) , 

(10) 
which is completely defined by the initial conditions and 
is independent of W{x,t). An initially purely right mov- 
ing wave packet, a„ = is again seen to generate an 
undistorted current density peak jx — v-p\a[x — vpi)\'^ /2 
moving at v^ towards the right. According to ([T0| . right 
and left movers in the initial wave will just add their con- 
tributions to the current density of opposite sign. Once 
more, this confirms the finding of perfect Klein tunneling 



through a barrier U{x,t) of any space and time depen- 
dence. Furthermore, the current jx does not depend on 
Tj, , giving the same contribution from both valleys K and 
K' . Regarding the current normal to the barrier wc find 
no unusual effects arising from the superposition of right 
and left moving amplitudes. 

Surprisingly, although we consider electron momenta 
with py = 0, we find a nonzero value for jy parallel to 
the barrier, 

jy{x,t) — TzVp\a+{x - Vpt)a^{x + Vpt)\ 

X Im(a+ale'(^+(^'*)-^-(^'*») . (11) 

Here, we define S± — Re(iS±) . This equation can be 
rewritten as 

jy(x,t) = r^Vp|a+a+(x,t)||a_a_(x,i) 

X sin[S+{x,t) - S^{x,t) + (j)] (12) 



with (j) = a,Tg{a+a*_). While eq. (fT2|) vanishes for uni- 
directional wave packets since then either a+ or «_ is 
zero, jy ^ for superpositions of left and right moving 
fermions. It is this superposition which causes the quali- 
tatively new phenomenon of a current along the barrier, 
exhibiting striking properties. The sine-dependence in 
(fT2|) is reminiscent of the Josephson effect where it orig- 
inates from the spatial overlap of superconducting order 
parameters in the contacts of a Josephson junction. In 
contrast to its a;-component, the y-component of the cur- 
rent manifests a nontrivial space and time dependence. 
The latter also depends on r^. Therefore, the best way 
to observe jy is to prepare a valley-polarized system. For 
non- or partly polarized situations one should add the 
current contribution from the second valley, yielding a 
total current jy{x, t) = jy{x, t)+jy {x, t) where jy and 
jy refer to states near K and K' points, respectively. 
Notice, that a contribution from the other Dirac point 
K' can compensate the current from K so that jy would 
allow to measure the degree of valley polarization in a 
graphene sample. However, even without valley polar- 
ization, the current fiuctuations in y-direction. 



(jy) {jy? = < 



1 



Ia+|2|a_|2 



(13) 



X SVC? {S+{x,t)- S-{x,t) + cj)) 



do not depend on t^, and therefore remain nonzero. 

Spatio-temporal dependence of jy . — To substantiate 
this effect we now investigate specific examples for po- 
tentials M^(a;,t) and the analogy between the current 
jy and a Josephson current for the valley-polarized sit- 
uation, Tz = 1. Let us consider the two potentials 
Wi{x) = WqxJL and W2{x) = WqCOs{xIV), both with 
characteristic length L and amplitude Wq. We consider 
periodic modulations in both cases with frequency a; so 
that W\2 = vanishes on time average. As initial con- 
dition we assume a superposition of right and left prop- 
agating plane waves, S'^ = ifcx . 




FIG. 1: (Color online) (a) Current jy (measured in units 
^•max _ ^p|q^q,_|^ along the barrier versus time, for the po- 
tential Wf'', see text, kx + 0/2 = 7r/8, and W^v^ j [Lb?) = 
I/tt (blue line) and W^v-pj^LiiP'') = 1/2 (red line), assuming a 
valley-polarized situation. For "Shapiro-step" conditions pe- 
riodic oscillations can be seen (red curve), while, away from 
this condition, aperiodic oscillations occur, (b) The same 
as (a) but for W{ potentials. This yields periodic solu- 
tions with a nonzero dc-current for any Wq {Wqv-pKLlJ''^ = 
I/tt (blue) or W^v-el(LuP') — 1/2 (red)) and no aperiodic 
regimes occur, (c) The same as in (a) but for potential 



W. 



(1) 



X — 'itL/2, k — 0, WqL/vf ~ 0.1, and frequencies 



Lj = {ti/2){vf/L) (red curve) and u = v^/L (blue curve). 
In both cases aperiodic oscillations occur. For the matching 
condition lo = v^/L (blue curve) a considerable enhancement 
followed by a saturation of the amplitude of the current oscil- 
lations occurs, while away of this resonance weak oscillations 
are not amplified. 



Substituting time-dependent potential IF^ {x,t) ~ 
Wo(a:/L) sinwt in eq. ^ and then using cq. P^ . we 
derive 



jy{x,t)=v^\a+\\a-\sm2 fcx + i^ (^ 



sin ujt ^ 



(14) 

for this particular potential. This equation can be rewrit- 
ten as a sum 



jy{x,t) 



WF|a+||a- 



f; ( 2WoVf \ 



n=0 
2WfiVFt 



X sin ( 2kx + Z1121l£. — |- nut + cj) ) (15) 

using Bessel functions J„. The last equation indicates 
that at uj — LUn with 



a;„ = y/2WaVF/{Ln) , neN 



(16) 



and integer n, the y-component of the current shows 
a peculiarity similar to the so-called Shapiro steps [28| 
of a Josephson junction. As seen in Fig. la, the fre- 
quencies uj = uJn generate periodic oscillations, which. 



again as in the case of Shapiro-steps, exhibit a nonzero 
dc-component in the current at given x. In view of 
eqs. (|14ll5p the overall dc-current vanishes after aver- 
aging over X. When modulating the potentials W^ ' 
with Lo ^ LOn results in aperiodic oscillations and zero 
rfc-component. 

The phase shifted modulation W^ (x,t) = 
Wq {x/L) cos ujt yields 

jyix,t) = wp|a+||a_| sin2 [fcx -f ^^W gi^s i^ + | 

(17) 
without a term proportional to t in the square bracket ar- 
gument of the sine- function, cf. eq. (|14l) . and thus without 
similarity to Shapiro steps in Josephson junctions. The 

(2) 

oscillations of jy generated by Wi are always periodic 
and with non-zero dc-component as seen in Fig. lb. The 
magnitude of this dc-current allows to directly measure 
the degree of valley polarization as needed for future val- 
leytronics. In the limit w ^ both versions for Wi 
become time-independent. 

The second choice, W2 ix,t) = Wq cos(a;/L) coswi , 
exhibits an even more intriguing time dependence due to 
possible spatio-temporal mode matching. In this case 



S+-S. 



2kx 



AWoLv^sinm 



cj^L^ — v'i 



(18) 



ujL - 



2L 



ljL - 



2L 



Now, even when w — > 0, the oscillations of jy persist 
since the static potential of spatial periodicity L induces 
a frequency component Vjr/ L to waves moving at the uni- 
form velocity v-p- This reminds of the ac- Josephson ef- 
fect 28] where ac-current oscillations are generated by a 

time-independent voltage. Shifting the phase of the time 

(2) 
dependence, W2 {x,t) — Wocos{x/L) sinujt produces 



S+ — S- — 2kx 



H^£|l!^(..Lsin(^)-..sinM)) 

(19) 
which vanishes when w — > 0. If spatio-temporal matching 
occurs, i.e. if a; — >■ vp/L, then both of the previous solu- 
tions vary proportional to t as 2kx =F t Wq sin ^ | ^J,"^* } 
for {g°^|^(} potentials. This results in an amplification 
of oscillations of jy with time for relatively small Wq, 
Fig. Ic. Indeed, even for small amplitudes of the scalar 
potentials, the oscillations grow resonantly at times t ~ 
2tt/Wo- We mention the analogy to the resonance excita- 
tion of plasmonic oscillations (Wood's anomaly [29J) by 
spatio-temporal matching of the incident light with the 
grating period. 

Concluding, we derived the exact solution of the Dirac 
equation for electrons in graphene moving strictly per- 
pendicular to a scalar barrier U{x,t), taken as homoge- 
neous along the j/-direction, which may vary arbitrarily 
with time and x. We still find unit transmission proba- 
bility, referred to as Klein tunneling, rendering at most 



a phase to the wave function. Importantly, we discover a 
nonzero current component jy parallel to the barrier, de- 
spite of vanishing incident momentum py . This current 
arises in valley polarized situations for packets containing 
both, left and right moving waves. It exhibits striking 
resemblance to the Josephson current of coupled super 
conductors and we have shown conditions under which 
Shapiro step like solutions arise. 
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